Optical properties of Ce02 using screened hybrid functional and GW-\-U methods 
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The optical spectra of Ce02 have been systematically investigated using three first-principles 
computational approaches for comparison, namely, the Heyd-Scuseria-Ernzerhof (HSE) screened 
hybrid functional, HSE+t/, and CW+U. Our results show that by using the HSE-I-/7 method, the 
calculated electronic structures are in good agreement with experimental spectra and the resulting 
imaginary part of the optical dielectric function spectrum well reproduces the main features exhibited 
in experimental observations. The important adsorption spectrum and energy loss function also 
accord well with the experimental results. 



Cerium oxides have been applied extensively because 
of the technologically important applications in industry 
originating from their unique physical properties. Taking 
Ce02 for example, it is used in automobile exhaust cata- 
lysts due to its capacity of high oxygen storage. Another 
novel application of Ce02 in the microelectronic and op- 
toelectronic fields has been attracted more attention be- 
cause of its high dielectric constant and good epitaxy on 
Si substrates resulted from the same cubic structure and 
small lattice mismatch Particularly, the dielectric 
function related optical properties of Ce02 have been in- 
vestigated intensively by several experimental work , 
and the complex 4/ orbitals of Ce are also proved to be 
unoccupied locating within the band gap formed between 
the O 2p and Ce 5d states. 

From theoretical viewpoint, accurate calculations of 
optical properties for 4/ rare-earth oxide systems are 
hard. For example, firstly, conventional density func- 
tional theory (DFT) within standard local or semilocal 
functionals always underestimates the band gaps of in- 
sulators or semiconductors, consequently, the resulting 
optical spectra are unsatisfactory. Secondly, an accurate 
description of electronic structures for rare-earth oxides 
containing 4/ electrons is also a great challenge to DFT 
due to the simultaneous itinerant and localized behav- 
iors exhibited by the / orbitals. Many efforts have been 
devoted to overcoming these intractable issues. 

One successful method to remedy the drawback of the 
standard local or semilocal functionals is the newly devel- 
oped hybrid functionals [5*1 which introduces 25% of the 
nonlocal Hartree-Fock (HF) exchange into the otherwise 
exact semilocal exchange functional. In order to apply 
this method to large molecules or extended systems, an 
improved version of the HSE screened Coulomb hybrid 
functional Q has been developed, in which the slowly de- 
caying long-ranged (LR) part of the HF exchange term is 
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replaced by the corresponding density functional counter- 
part. Taking Perdew-Burke-Ernzerhof (PBE) potentials 
for example, the resulting exchange is given by 

where E^^'^^ is the short-ranged (SR) HF exchange, 
gPBE,SR aj^j gPBE,LR ^j,g ^j^g gj^ ^.^^ components 

of the PBE exchange functional, respectively. 

Another successful method for determination of ex- 
cited states is the GW approximation |7|, which is based 
on the quasiparticle concept and the Green function 
method. In this approximation, the self-energy S is ex- 
panded linearly in terms of the screened interaction W, 
i.e., S w GW. Here, G is the single-particle Green's 
function and W is the screened Coulomb interaction. In 
order to describe the localized / orbitals more accurately, 
we also adopt the PBE plus a Hubbard U correction 
(PBE-I- U) method as the starting point, which gives 
a qualitative improvement compared to PBE not only for 
excited states such as band gap but also for ground-state 
properties such as Mott insulator. 

In this letter, we systematically investigate the optical 
spectra of CeOz using the HSE, HSE-I-C/ and GW+U 
methods as implemented in the first-principles vienna 
ab initio simulation package (VASP) [9|. The electron 
and core interactions are included using the frozen-core 
projected augmented wave (PAW) approach ,10,] . The 
gradient corrected PBE functional [ll| for the exchange 
correlation potential is used. The cerium 5s, 6s, 5p, 5d 
and 4/ as well as the oxygen 2 s and 2p electrons are ex- 
plicitly treated as valence electrons. The electron wave 
function is expanded in plane waves up to a cutoff en- 
ergy of 500 eV. For the Brillouin zone integration, the 6 
X 6 X 6 F-centered fc-mesh is adopted and a good con- 
vergence can be obtained. In the HSE calculations, 25% 
HF exchange and 75 % PBE exchange are used, and the 
range-separation parameter is set to ^=0.3 A~^. The 
GW calculations are carried out using a total number of 
166 bands and four iterations for accurate quasiparticle 
shifts. The strong on-site Coulomb repulsion among the 
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localized Ce 4/ electrons is described by the rotationally 
invariant approach to the LDA + U method due to Du- 
darev et al. [l^. In our calculation, we use J — 0.51 eV 
for the exchange energy and U — 4.5 eV for the spheri- 
cally averaged screened Coulomb energy which is close 
to the values used in other previous work [l3| . Since only 
the difference between the spherically averaged screened 
Coulomb energy U and the exchange energy / is im- 
portant for the total LDA (GGA) energy functional [l^ . 
thus, in the following we label them as one single effective 
parameter U for brevity. Notice that all our calculations 
are performed at the experimental lattice constants. 

Before calculating the optical spectra of Ce02 , we first 
study the essential electronic structure using the HSE, 
HSE-f {/ and GW+U methods, respectively, since the 
spectra are calculated directly from the interband tran- 
sitions. Based on this, the frequency-dependent complex 
dielectric function e{uj) — ei{uj)+ie2{uj) for Ce02 are cal- 
culated. In the following, we will show and discuss our 
results systematically. The calculated results obtained 
using standard DFT are not present in this letter, be- 
cause they can be available elsewhere p"3l - [l5| . Figure 
1(a) shows the site-projected orbital-resolved density of 
states (DOS) for Ce 4/, 5d and O 2p orbitals using HSE 
method together with the experimental XPS + BIS spec- 
tra [IBl- From the DOS, it is clear that the valence bands 
mainly consisted of O 2p states and a little contribution 
from Ce 4/ and 5d states can also be seen. The con- 
duction bands mainly consisted of Ce 5d states, and the 
localized Ce 4/ states are located within the p—d gap. 
Our calculated width of O 2p valence bands is about 
4 eV, which is in good agreement with the experimen- 
tal XPS results [la]. Note that the fundamental p—f 
and p—d gaps govern the optical transitions, therefore, 
their accurate determination is indispensable to the op- 
tical spectra. Our calculated p—f gap is about 2.8 eV, 
which is well consistent with the experimental value of 
3.0 eV [ll]. Furthermore, the p—d gap of about 6.2 eV 
is also well determined compared with the experimental 
value of about 6.0 eV [l^. Another remarkable property 
is the central position of localized Ce 4/ orbitals. Our 
DOS showed in Fig. 1(a) shows that the localized Ce 4/ 
orbitals lie about 3.5 eV above the topmost of the O 2p 
valence bands, which is smaller than the XPS-BIS result 
of about 4.5 eV [lB|. Notice the that present electronic 
structure obtained by HSE arc in good agreement with 
previous study |l5| using HSE. 

Based on the determined electronic structure, the op- 
tical spectra of the frequency-dependent dynamical di- 
electric function s{uj) = £i(w) -I- K2(i^) as a function of 
the photon energy uj up to 12 eV has been calculated 
and showed in Fig. 1(b). In order to make a compari- 
son, the experimental result [2] of imaginary part £2(1^) 
is also showed. Due to the cubic symmetry of Ce02, 
the dielectric tensor only has one independent compo- 
nent and Sxx = £yy = £z2 7 therefore, only one component 
is showed in Fig. 1(b). According to our above DOS dis- 
cussion, the peak at about 3.9 eV in our calculated £2('^) 



is assigned to the transition from the O 2p bands to the 
localized Ce 4/ states, which is close to the peak in ex- 
perimental £2(w) at 3.7 eV [1]. Considering the dipolar 
selection rule only transitions with the difference = 
±1 between the angular momentum quantum numbers I 
are allowed, the p—f transition is forbidden. However, 
in our above DOS calculation, in valence bands, small 
hybridization between O 2p and Ce 4/ exists, which sug- 
gests that some Ce 4/ electrons transfer to O 2p states. 
Notice that Niwano et al. Q also called the p—f tran- 
sition as a "charge transfer" transition. As for the O 
2p—Ce 5d transition, experimentally Marabelli et al. Q 
observed a double-peaked very broad structure at about 
9 eV in £2(w); one peak is centered at 8 eV and the other 
at 10 eV. This corresponds to the transitions from O 2p 
to Ce 6d eg and t2 states due to crystal field splitting in 
fluorite crystal structure. As for the p—d transition indi- 
cated in £2(w) in Fig. 1(b), one can see two peaks located 
at about 8.6 and 10.8 eV, corresponding to p—d{eg) and 
p—d(t2) transitions, respectively. 

The electronic structure and optical spectra obtained 
by the HSE -|- U method are showed in Fig. 1(c) and Fig. 
1(d). Compared to our pure screened hybrid functional 
result, it is clear that the Ce 4/ states are more localized 
because of the Hubbard U correction. Furthermore, the 
central position of localized Ce 4/ states lies at about 4.5 
eV above the topmost of the O 2p valence bands, which 
becomes to accord well with the experimental XPS-BIS 
measurement [1^1 . As for the p—f and p—d gap s, they 
are now very close to the experimental data |l6]. Over- 
all, the shape of our calculated DOS curve exhibits the 
same main features as the experimental XPS-BIS results 
[l6| . Correspondingly, the calculated £2(0;) by the HSE 
-I- [/ is showed in Fig. 1(d). Clearly, there are three 
peaks in the calculated £2(1^)- One peak is at 4.9 eV, 
and the other two are at 8.1 and 10.3 eV, respectively. 
Comparing to the three main peaks (located at 3.7, 8 
and 10 eV) observed by experiment we notice that 
the last two peaks are well reproduced except the first 
peak with a shifting of about 1.2 eV from the experi- 
mental spectrum. This discrepancy is guessed to be due 
to the ignorance of the interaction between electron and 
hole in our calculations. Then, the excitonic effect is 
not taken into account. As Niwano et al. Q pointed 
out that the BIS process is adding one additional elec- 
tron into the empty conduction bands, whereas optical 
adsorption process is exciting one electron up to the con- 
duction bands and a corresponding hole is created; they 
further deduced the binding energy of an exciton con- 
tributed by one hole and one 4/ electrons to be about 
1 eV according to energy difference between the central 
position of the localized 4/ states and the first peak in 
£2. Notice that compared to the experimental observa- 
tion [1], our £2 spectrum also exhibits a rising edge at 
about 6.5 eV, which is close to the p—d gap and indi- 
cates the beginning of the p—d transitions. Therefore, 
our calculated £2 spectrum displays the main features as 
the experimental result and supplies a qualitative agree- 
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FIG. 1: The DOS and dielectric function spectra along with corresponding experimental results for Ce02: (a), (c), and (e) the 
projected orbital-resolved partial DOS for O 2p, Ce 4/ and Ce 5ci orbitals using HSE, HSE + U, and GW + U, respectively; 
(b), (d), (f) the optical spectra of the frequency-dependent dynamical dielectric function e{(jj) = ei(i^) + K2(i^) as a function of 
the photon energy u using HSE, HSE -I- U, and GW + U, respectively. The experimental XPS-BIS and £2 spectra are taken 
from Refs. J£] and j2], respectively. 



mcnt with experimental spectra [2, l3|. Furthermore, if 
the excitonic effect is taken into account, the oscillator 
strength in 62 will shift to lower energies, we expect that 
the agreement will become better. 

The calculated projected DOS and dielectric function 
spectra obtained by the GW+ U method are also showed 
in Fig. 1(e) and Fig. 1(f). It is clear that the Ce 4/ 
orbitals are further localized compared to the HSE+ U 
results. As for the gaps, the p—d gap is close to exper- 
imental result, whereas the p—f gap of about 3.9 eV is 
enlarged. Notice that the band gaps for semiconductors 
or insulators are usually overestimated in GW scheme. 
This is because the screening is reduced upon updating 
the eigenvalues in W [13]. As expected, the first peak 
in £2 spectrum will shift to a higher energy and our cal- 
culated £2 spectrum shows that it is located at 6.0 eV, 
while other two peaks at 8.1 and 10.3 eV, respectively. 

Comparing the three results showed in Fig. 1, one 
can conclude that the results obtained by the IISE+ U 
method are better than the other two. Furthermore, for 
£2 spectrum, we expect the agreement will be become 
better if the excitonic effect is taken into account. Based 
on this, the important adsorption spectrum and energy- 
loss function deduced from the frequency-dependent dy- 
namical dielectric function obtained by HSE-I- U, together 
with the corresponding experimental data Q , are showed 
in Fig. 2. It is obvious that the main structures are well 
reproduced. Similar to the £2 spectrum, the first peak 
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FIG. 2: The HSE-I- U results together with corresponding ex- 
perimental results for Ce02, (a) the optical adsorption spec- 
trum, (b) the energy loss function: orange solid line is this 
work, blue short dash dot line is experimental results taken 
from Ref. jl]. 



in adsorption spectrum has a shifting of about 1.2 eV to 
higher energy. As for the energy-loss function, the two 
peaks at about 5.5 and 15.5 eV also shift to high energies 
by about 1.1 eV compared to the experimental results 
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In summary, we have studied the optical properties of 
CeOa using HSE, HSE+f/, and GW+U three methods. 
Using HSE+ U method, our calculated band gaps and 
the central position of 4/ orbitals are in good agreement 
with the experimental values, and the resulting optical 
dielectric function, adsorption spectrum and energy-loss 
function obtained by HSE-|- U are also well reproduced. 
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